Mon. Not. R. Astron. Soc. 0Q0.[T1-fTT1(20051 Printed 2 February 2008 (MN teTpK style file v2.2) 



Seismic analysis of the second ionization region of helium in the 
Sun: I. Sensitivity study and methodology 

Mario J. P. F. G. Monteiro 1 *' 2 and Michael J. Thompson 3 * 

1 Centra de Astrofi'sica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

2 Departamento de Matemdtica Aplicada, Faculdade de Ciencias da Universidade do Porto, Rua do Campo Alegre 687, 4169-007 Porto, Portugal 

3 Department of Applied Mathematics, University of Sheffield, Sheffield S3 7RH, UK 



Accepted 2005; Received 2005; in original form 1996 



ABSTRACT 

The region of the second ionization of helium in the Sun is a narrow layer near the sur- 
face. Ionization induces a local change of the adiabatic exponent T-y, which produces a char- 
acteristic sign ature in the frequencies of p-modes. B y adapting the method developed by 
Monteiro, Christensen-Dalsg aard & Thompson! i 1994ft . we propose a methodology for deter- 
mining the properties of this region by studying such a signature in the frequencies of oscilla- 
tion. 

Using solar data we illustrate how the signal from the helium ionization zone can be 
isolated. Using solar models with different physics - theory of convection, equation of state 
and low temperature opacities - we establish how the characteristics of the signal depend on 
the different aspects contributing to the structure in the ionization layer. We further discuss 
how the method can be used to measure the solar helium abundance in the envelope and to 
constrain the physics affecting this region of the Sun. 

The potential usefulness of the method we propose is shown. It may complement other 
inversion methods developed to study the solar structure and to determine the envelope helium 
abundance. 

Key words: Sun: interior - Sun: helioseismology - Sun: oscillations - Sun: abundances - 
equation of state - stars: abundances 



1 INTRODUCTION 

The direct determination of the helium abundance in the solar near- 
surface layers is difficult and uncertain, although it is very impor- 
tant to the modelling of the internal structure and evolution of the 
Sun ("see lKosovichev et alll 992 for a comprehensive discussion of 
the sources of uncertainties). It is usually taken as a fitting parame- 
ter of an evolutionary sequence that provides the correct luminosity 
for the Sun at the present age. The possibility of constraining this 
parameter to have the observed value for the Sun is important to 
improve the mass loss estimates and early evolution of the Sun, as 
well as to test the effects of mixing and settling on stellar evolution. 

Several attempts have been made to use solar seis- 
mic data to calculate the abundance of helium (Y) in the 
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poses serious difficulties to an accurate direct seism ic measure- 
ment of t he envelope abundance of helium jKoso vichev et ail 

^J92| [PerezHernandez & Qtristensen-Daisgaa^ "" Hl^T 

Basu & Christensen-Dalsgaard 1997). The sensitivity of the 
modes to the helium abundance is primarily provided by the 
change of the local adiabatic sound speed c due to ionization. Such 
sensitivity is given by the behaviour of the first adiabatic exponent, 
Ti, since c 2 =Fip/p where p and p are the pressure and density 
respectively, and consequently it is strongly dependent on the 
assumed equation of state and other physics relevant for the region 
where the ionization takes place. This is the main reason why the 
seismic determination of the envelope abundance of helium is 
highly complex. 

Here we propose a method complementary to those 
use d previously, by a dapting the procedure developed 
bv iMonteiro et alJ Jl994l in the follow in g M CDT) and 
Christensen-Dalsgaard, Monteiro & Thompson 1 1995, in the 
following CDMT). In using the solar frequencies in a different 
way, which provides a direct probe to the region of ionization, we 
aim at providing a method where the different effects at play in the 
ionization zone can be isolated, constructing a procedure to access 
the chemical abundance. Localized variations in the structure of 
the Sun, such as occur at the base of the convective envelope 
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(see MCDT and lMonteirclll99d) and in the region of the second 
ionization of helium lMonteiroll996h . create a characteristic signal 
in the frequencies of oscillation. The properties of such a signal, as 
measured from the observed frequencies, are related to the location 
and thermodynamic properties of the Sun at the layer where the 
sharp or localized variation occurs. The main advantage we see 
in this method is the possibility to utilise different characteristics 
of the signal to distinguish different aspects of the physics of the 
plasma at the region where the signal is generated. In particular 
we may be able to separate the effects due to convection, the 
low-temperature opacities and the equation of state from the quan- 
tification of the helium abundance that we seek to achieve. Here 
we mainly concentrate on separating these different contributions 
in order to establish the dependence of the parameters of the signal 
in the frequencies on the different aspects of the structure at the 
ionization region. Using a variational principle we determined 
how the zone of the second ionization of helium can indeed be 
considered as a localized perturbation to an otherwise 'smooth' 
structure, generating a characteristic signal in the frequencies of 
the modes. 

We note that simplified versions of the expression for the 
signal discussed here have been applied successfully to cases 
where there are only very low degree frequencies. The signal 
has been fitted either to the frequencies of low degree modes 
iMonteiro & Thompsoij|l998tlVerner. Chaplin & ElswortrJl2004h. 
or to frequency differenc es (Miglioetal. 2003; Basu et al. 
1 2004 IVauclair & Theadol l20ol iBazot & Vauclairl 120041 
IPiau. Ballot & Turck-Chieza l2005t) . Here we obtain the ex- 
pression for the general case of having also modes of higher 
degree, of which the low degree applications are a particular case. 
We also demonstrate the method for making use of the information 
in moderate-degree data available only for the Sun. When using 
modes with degree above 4 or 5 we can avoid using frequencies 
affected by the base of the convection zone and may hope to 
achieve a much higher precision in the results as many more 
frequencies with lower uncertainties can be used. 

In this work we present the analysis of the characteristics of 
the signal under different conditions. Several models with differ- 
ent physics and envelope helium abundances are used to test the 
method in order to prepare the application to the observed solar 
data. 



2 THE REGION OF THE SECOND IONIZATION OF 
HELIUM 

In order to model the sensitivity of the modes to this region we 
must try first to understand how ionization changes the structure. 
Secondly, we need to estimate how the modes are affected by such 
a region. The details of the derivations are discussed in the Ap- 
pendix, but the assumptions and the main expressions are reviewed 
and analysed here. 

2.1 Properties of the ionization region 

Because the helium second ionization zone (Hell ionization zone) 
is sufficiently deep (well within the oscillatory region of most of 
the modes) we propose to adapt the method discussed in MCDT 
to the study of this layer. The contribution from a sharp variation 
in the structure of the Sun to the frequencies can be estimated by 
calculating from a variational principle for the modes the effect of 
a localized feature. In the work by MCDT the feature was the base 
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Figure 1. Plot of the adiabatic exponent Ti for several solar models. As 
a reference we have calculated a model (Zq) where the second ionization 
of helium has been suppressed. The other two models are calculated using 
different equations of state (see Table 1 for further details on the models). 
The second ionization of helium takes place around an acoustic depth of 
600 s, corresponding to the depression on the value of Ti . 

of the convection zone and the sharp transition was represented by 
discontinuities in the derivatives of the sound speed. Because of the 
size of the ionization region when compared with the local wave- 
length of the modes, that representation is not adequate to repro- 
duce the effect on the frequencies for the ionization region. 

Here we must, instead, consider how the ionization changes 
the first adiabatic exponent ri=(d\np/dln p) s (the derivative at 
constant specific entropy s) locally, generating what can be de- 
scribed as a 'bump' over a region of acoustic thickness of about 
300 s (see Fig.0. This allows us to estimate how the frequencies 
of oscillation are 'changed' due to the presence of this feature in 
the structure of the Sun. The effect will be mainly taken into ac- 
count through the changes induced in the adiabatic gradient Ti by 
the ionization. Other thermodynamic quantities are also affected, 
but the changes on the local sound speed is mainly determined by 
changes in Ti. Therefore, we will calculate the dominant contribu- 
tion to the changes in the frequencies by establishing what is the 
effect on the modes due to changes of the adiabatic exponent. 

| E>appen & Goughl ll986l) and iDappen. Gough & Thompsor] 
(1988) have proposed a method based on the same principle, by 
using the sensitivity of the sound speed to changes on the adia- 
batic exponent. Using this sensitivity they calibrate a quantity that 
is associated with ionization in order to try to measure the helium 
abundance in the solar envelope from seismic data. But most meth- 
ods have difficulties in removing the dependence of the calibration 
on the physics of the reference models, making it difficult to obtain 
a value for the abundance. This is the problem we try to address 
in this contribution by developing a method able to measure in the 
frequencies the effect of the ionization and its dependence on the 
abundance, isolated as much as possible from the other uncertain- 
ties. 

2.2 A variational principle for the effect on the frequencies 

A variational principle for nonradial adiabatic oscillations, assum- 
ing zero pressure at the surface located at radius R as a bound- 
ary condition, can be formulated. It is possible to further consider 
only higher-order acoustic modes, for which we may neglect the 
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Figure 2. (a) Plot of the differences (<5ri/ri) between two models - one 
with and the other without the second ionization of helium, versus acous- 
tic depth t. These correspond to models Zq and Z\ discussed in the text 
and described in Table 1. The dotted line represents our assumed smooth 
reference structure, (b) The change of V\, relative to the smooth reference 
structure, is shown. The parameters represented schematically, to be deter- 
mined from the characteristics of the signal in the frequencies, are: 8^ - the 
amplitude of the differences at r^, and /3 - the half width of <5ri/ri (values 
taken from Tablel3l below). Also indicated is the value of f c i as found from 
the frequencies. 



perturbation in the gravitational potential. The outcome of such a 
variational principle is an estimate of how the frequencies change 
due to changes in (Tip) and p. This is described and discussed in 
Appendix IaI 

In order to model the signature of the ionization zone we rep- 
resent the effect of the second ionization in terms of the changes it 
induces in the adiabatic exponent Ti . Such a change (see Fig|2j is 
approximately represented by a 'bump' of half width f3 in acoustic 
depth, and relative height 
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being located at a radial position corresponding to an acoustic depth 
Td. Here, and in the following, acoustic depth r at a radius r is 
defined as, 



r(r) 



dr 

c 



(2) 



where R is the photospheric radius of the Sun. 

Relatively to the frequencies of a reference model, assumed to 
be 'smooth' and corresponding approximately to a model with no 
Hell ionization region, we find that the bump changes the frequen- 
cies in such a way that there is a periodic component of the form 
(see AppendixlAl. 



Suj ~ A(uj, I) cos Ad 



(3) 



where the amplitude, as a function of mode frequency lo and mode 
degree I, is given by 



A(ui, I) = ao 



1-2A/3 sin 2 [/3w(l-A) 



1/2] 



(1-A) 2 
and the argument is 



/3w 



A d 



(1-A) 1/2 dr + 



(4) 



(5) 



Here the factor in A represents the geometry of the ray-path ac- 
counting for deviation from the vertical when the mode degree is 
non-zero. It is associated with the Lamb frequency, as given below 
(Eqs[8]|9j. In fact, because the ionization zone is close to the sur- 
face and provided we are not using very high-degree data, we can 
neglect A in the expression for the argument Ad; we can similarly 
neglect the effect of the mode degree on the surface phase func- 
tion <f>. Consequently, for the ionization zone the expression of the 
argument becames 



A d - 2 (W d + <f>) ~ 2 (ur d + 0o ) 



(6) 



In the asymptotic expression used for the eigenfunction (see 
Eq. lA4t . the phase 4> depends on the mode frequency and degree 
(see MCDT for details). Here we have expanded the phase to first 
order in frequency, by writing that (j!>(w)~</>o-|-a0tj. From this it 
follows that rd=Td+a^, while the frequency independent term of 
the phase is now 4>o- 

The amplitude of the signal, through ao, corresponds to 
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where r t =r (0) is the total acoustic size of the Sun. The small factor 
A, present in the amplitude, is given by 



A = A c 



1(1+1) Co 2 



1(1+1) W 
where the value of Ad is given by, 
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and I and u are two reference values. These values are chosen tak- 
ing into account the expected behaviour of the signal and the set of 
modes used, as discussed below. 

In order to compare the amplitude it is convenient to define a 
reference value Ad, as given by 



Aa 



A(Cj,1) 
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ao 
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0Cj 



(10) 



The parameters of the signal relevant to characterize the local 
properties of the ionization zone, as given in Eq. are f<j, (3, ao 
and Ad- 

The values of fd and Ad can be used to measure mainly the 
location of the ionization zone. They both vary strongly with dis- 
tance to the surface. But the acoustic depth is a cumulative function 
of the sound speed behaviour over all layers from the surface to a 
particular position, whereas Ad is a local quantity, not being af- 
fected by the behaviour of the sound speed in the layers above the 
ionzation zone. 

The values of and ao (or Sd) are expected to be directly 
related to the local helium abundance, since the size of the bump 
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will be determined by the amount of helium available to be ionized. 
These parameters are also expected to be strongly affected by the 
equation of state, and to a lesser extent by the other physics that 
affect the location of the ionization zone (ra). But we may hope 
to be able to use the other parameters to remove this dependence, 
while retaining the strong relation between the bump and the he- 
lium abundance (Y). 

2.3 Measuring the signal in the frequencies 

Our first goal is to find the five parameters describing the signal 
from the frequencies of oscillation. In order to do that we must 
isolate a signature of about 1 /iHz in amplitude, overimposed in 
actual frequencies. We do so by isolating in the frequencies the 
periodic signal described by Eq. |3} using a non-linear least-squares 
iterative fit to find the best set of parameters. The method used is an 
adaptation of the one proposed by MCDT; for the present problem 
we must redefine the characteristic wavelength of the signal to be 
isolated (quantity Ao in MCDT) as it is significantly larger than for 
the signal from the base of the convective envelope. The parameters 
describing the signal (Eq.|3J, and found by our fitting procedure, are 
the following; 

Ta, 0o, do, Ad, (3 . 

We choose a set of modes which cross the ionization zone, but 
which do not cross the base of the convection zone. By remov- 
ing modes that penetrate deep in the Sun (low degree modes), we 
avoid the contamination coming from the signal generated at the 
base of the convection zone (see MCDT). But when selecting only 
modes of higher degree (between 45 and 100), it becomes nec- 
essary to include the contribution from the mode degree to the 
amplitude of the signal. This is the reason why it is necessary 
to include in the fitting the paramet er Ad. This parameter is not 
necessary when studying other stars l Monteiro & Thompson 1998; 
iBasu et all20ol ; |Piau et alj20 05). resulting in a simplified descrip- 
tion of the expected observed behaviour. In the case of the Sun there 
is a great advantage in using all available high-degree modes that 
cross the ionization zone. 

The modes considered correspond to the ones available in so- 
lar data, having degrees and frequencies such that the lower turning 
point is between 0.75R and 0.95R of the solar radius. The latter en- 
sures the modes cross the ionization zone while the former avoids 
contamination from the signal originating at the base of the con- 
vective envelope (e.g. CDMT, and references therein). These con- 
ditions define typically a set of about 450 modes having frequency 
lu/2h in the range [1500, 3700] /iHz, and with mode degree of 
45^100. 

As we are only using modes of high degree in this work, the 
reference values preferred in the fitting of the signal are; 

I = 100 and — = 2000 ttHz . 

2tt p 

The first value is an upper limit for modes that cross beyond the 
ionization zone, while the value of Co corresponds to the region in 
frequency where the signal is better defined. These values are only 
relevant to normalize the parameters fitted for different models. 

For solar observations only frequencies with a quoted observa- 
tional error below 0.5 /iHz are included. We ensure consistency of 
the data sets by restricting the selection of mode frequencies from 
the models to the modes present in the solar data after applying the 
above selection rules. 

We stress that the method adopted for removing the smooth 
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Figure 3. (a) Residuals left after a smooth component of the frequencies, 
as functions of mode order and degree, is removed. The data are from 
MDI/SoHO, with the error bars corresponding to a 3a of the quoted ob- 
servational eiTors. (b) Plot of the signal isolated, and shown in panel (a), 
after division by the amplitude function as given by Eq. J4} when using the 
values of (ao, <5di /?, Ad) found in the fitting. Modes with degree below 
£=60 are shown as filled circles, while modes with a higher value of the 
degree are represented by open circles. The line indicates the fitted periodic 
signal as expected from Eq. 

component of the frequencies is a key assumption in the process 
of fitting the signal. Here we use a polynomial fit with a smooth- 
ing parameter on the third derivative (see CDMT). In any case, as 
long as the analyses for different models and for the solar data are 
consistent, the comparison of the parameters will be independent 
of the choice on how to describe the smooth component. Such con- 
sistency is ensured by using exactly the same set of frequencies and 
the same numerical parameters for the fitting. 

2.4 The signal in the solar data 

To illustrate the signal extraction, the method proposed here was 
applied to the ana lysis of solar seism ic data from MDI on the 
SOHO spacecraft JScherrer et aflll995l) . The signal was isolated 
as described above for the models. The fitted signal of the Sun is 
shown in Fig. [3k together with the error bars. In order to illustrate 
how well the expression for the signal (Eq. |3J fits the data points 
we also show in Fig.|5J) the signal in the frequencies normalized by 
the amplitude as given in Eq. J4j. The quality of the fit done with 
Eq. confirms the adequacy of the first order analysis developed 
in Appendix[A]leading to the expression given by Eq. J4j. 

The values of the parameters found for the data are given in 
TableQ From Monte Carlo simulations we have estimated the un- 
certainty in the determination of the parameters due to observa- 
tional uncertainties as indicated by the quoted observational errors. 
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Table 1. Parameters obtained by fitting the observed solar frequency data 
with the expression of the signal as given in Eq. |3j- The quantities and f} 
are given in seconds, while the amplitudes (ao and A d ) are given in fiHz. 
Note that is not a fitting parameter, as it is given from the other pa- 
rameters using Eq. 1101 . The standard deviations cr are estimated from 200 
simulations of the effect of the observational uncertainties on the determi- 
nation of the parameters. 
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The values found, at the 3a level, are also listed in Table Q It is 
clear that due to the large amplitude of this signature (above 1 pHz 
at w/2-7r=2000 /iHz) the precision with which the parameters are 
determined is very high. As long as the method to isolate this char- 
acteristic signature is able to remove the "smooth" component, the 
results can be used with great advantage for testing the zone of the 
second ionization of helium in the Sun. 



3 SOLAR MODELS WITH DIFFERENT PHYSICS 

In order to establish how sensitive the different characteristics of 
the signal are to the properties of the ionization zone, and there- 
fore to the aspects of the Sun which affect the ionization zone, we 
consider different static models of the Sun calculated with the same 
mass, photospheric radius and luminosity. The profile of the helium 
abundance in the models is obtained by calibrating with a constant 
factor a prescribed abundance profile from an evolution model with 
the age of the Sun (without settling). 

We note that imposing the same radius and luminosity for all 
models is the key difference between the analysis presented here 
and the work bv lBasu et alj 120041) . If the models are not required 
to have the same luminosity and radius as the Sun, the properties of 
the ionization zone are not affected in the same way. Consequently 
the behaviour of the amplitude of the signal in this case is different 
from what we find when both these conditions are imposed on the 
models. 

The aspects of the physics being tested here are the equation 
of state (EoS), the theory of convection and the opacity. All these 
aspects affect the ionization zone by changing its location, size and 
thermodynamic properties. 

All models wer e calculated as in Monteiro ll996t see also 
iMonteiro et alj [l996V These are not intended to represent accu- 
rately the Sun, but simply to illustrate the usefulness of the method 
we propose to study a particular region of the solar envelope. 

As the simplest possible EoS we have used a Saha equation 
of state with full ionization at high pressure - this corresponds to 
SEoS in Tabled As a more complete EoS we have used the CEFF 
equation of state as described in Christensen-Dalsgaard & Daooen 
( 1992). For the opacities we have considered a simple power law 
fit (SOp), or the Rosseland mean opacity tables at low temperatures 
from lKuruczNl99ll) . To include convection we have taken the stan- 
dard mi xing length theory iBohm-Vitensd Il958t parametrized as in 
Monte iro et aljl9 96) or the more recent CGM model ICanuto et alJ 
Il996t 

As our reference model, in order to illustrate the changes due 
to the ionization of helium, we have calculated a very simple solar 
model (Zo) with suppressed Hell ionization, by setting the ioniza- 



Table 2. Solar models and their helium (Y) abundances. Also indicated are 
the equation of state (EoS): SEoS - Simple Saha equation of state with pres- 
sure ionization, and CEFF (cf. Christensen-Dalssaard & Dannen 1992); the 
Opacity: SOp - simple power law fit o f the opacities, and Kur - low tem- 
perature opacities from Kurucz 1 1991); and the formulation for modelling 
convection: MLT - standa rd mixing length theory |Bohm-VitenseHl958J; 
parametrized as in lMonteiro et alll996l) . and CGM - ICanuto etal. (1996). 
See the text for a description of the parameter f t used in the calculation of 
the emissivity. 
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Figure 4. Plot of the differences for Ti between all models considered and 
the one without the second ionization of helium. See Table l2l for the de- 
tails of each model. Only the region around the second ionization of helium 
is shown corresponding to the negative bump around an acoustic depth of 
650s. 



tion potential to zero. The helium abundance found for each model 
corresponds to the value that fits the boundary conditions (by scal- 
ing a prescribed dependence of the chemical profile, as taken from 
an evolved solar model). 

The behaviour of the adiabatic exponent for some of the mod- 
els (see Table|2j, relative to our reference model (Zo), is illustrated 
in Fig.|4] There is a clear difference on the location of the ionization 
zone (r d ) when a different EoS is used. The effects of changes in 
the formulation of convection or in the opacities are much smaller. 

In order to have models with the same envelope physics, but 
different helium abundances, we have calculated solar models with 
the energy generation rate changed by a prescribed factor / e in the 
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Figure 5. (a) Plot of the adiabatic exponent Ti for all models '5', calcu- 
lated with the same physics but different surface helium abundances as im- 
posed by / e (see Tablel2*l. The hydrogen and both helium ionization regions 
are shown, (b) Plot of the differences for Ti between models Z5 5^ 5,, 
and model Z 5 ;. See Tablel2lfor the details of each model. Only the region 
around the second ionization of helium (indicated by the arrow) is shown. 



emissivity. These are models Zzi^h and Z^^h^v which are sim- 
ilar to Z3 and Z5, respectively, except for the value of / c which 
is now different from unity. These correspond to models with a 
different structure of the core but with envelopes with exactly the 
same set of physics. All differences between these models in the 
envelope are due to differences in the chemical composition. To 
illustrate the differences we plot in Fig.[3]the differences in Fi be- 
tween models with the same physics but increasing values for the 
envelope abundance of helium. As the helium abundance increases, 
there is a corresponding decrease in hydrogen, which results on a 
slight separation in temperature of the three major ionization re- 
gions. Consequently both ionization regions for the helium expand 
towards higher temperatures. As the bump becomes slightly wider 
and moves to a higher temperature, the effect on the frequencies is 
expected to become smaller. 

For all models we have calculated the frequencies of linear 
adiabatic oscillations. The set of frequencies for each model, as 
used to fit the signature of the ionization zone, is described above. 
The parameters obtained in fitting Eq. |3} to the frequencies of the 
models listed in Table[2](excluding Zq) are given in Table|3| 



Table 3. Parameters obtained by fitting the frequency data for the models 
with the expression of the signal as given in Eq. |3j. The quantities Td and 
f) are given in seconds, while the amplitudes (ao and A<j) are given in 
/iHz. Note that is not a fitting parameter, as it is given from the other 
parameters using Eq. 1101 . 



Model 




00 


ao /2n 


A d /2n 


P 


A d 


Zi 


718.0 


2.588 


1.634 


2.834 


142.5 


0.604 


z 2 


724.8 


2.525 


1.671 


2.862 


141.9 


0.599 


Z31 


729.9 


1.950 


2.500 


3.251 


146.0 


0.484 


Z3 


730.4 


1.951 


2.380 


3.140 


144.7 


0.490 


Zzh 


730.4 


1.951 


2.314 


3.066 


144.3 


0.491 


z± 


739.9 


1.859 


2.353 


3.151 


143.3 


0.495 


Z51 


737.7 


1.874 


2.429 


3.241 


143.7 


0.494 


Zs 


737.8 


1.876 


2.342 


3.145 


143.1 


0.496 


Z$h 


737.5 


1.880 


2.278 


3.072 


142.7 


0.498 


Zg v 


736.8 


1.890 


2.205 


3.002 


141.7 


0.502 


z (i 


746.4 


1.790 


2.280 


3.141 


141.3 


0.507 



4 THE EFFECT OF THE PHYSICS ON THE 
CHARACTERISTICS OF THE SIGNAL 

The set of solar models considered here, and listed in Table|2| cover 
three major aspects of the physics which determine the surface 
structure of the models: the equation of state, the low temperature 
opacities and the formulation for convection (defining the supera- 
diabatic layer). In order to use the diagnostic potential of this char- 
acteristic signature in the frequencies we need to understand how 
each parameter describing the signal is affected by a specific aspect 
of the physics defining the structure of the envelope. 

One would expect that the shape of the bump is strongly de- 
termined by the EoS. But any change in the structure that affects 
the location of the ionization zone will necessarily introduce an ef- 
fect on the characteristics of the Ti profile. Consequently we need 
first to identify the parameters that depend more strongly on the 
location. These are most likely fa and Ad- 

The changes on the upper structure of the envelope are ex- 
pected to have a direct effect on the turning point of the modes. 
Consequently we need to look at the parameters that may be af- 
fected by the upper reflecting boundary. This is mainly expected to 
be 0rj. 

Finally the area of the bump in Ti in the ionization zone 
should reflect the local abundance of helium, if the location is well 
defined. Therefore we will look at ao and (3 in order to identify how 
the helium abundance Y defines the characteristics of the signal in 
the frequencies. 

4.1 The location of the ionization zone 

The most easily identifiable characteristic of the signal is its pe- 
riod. This quantity depends strongly on Td, but as discussed when 
writing Eq. J5J the period also contains a contribution from the up- 
per turning point of the modes (where there is a phase shift of the 
eigenfunction). This means that the period, or precisely fd, that we 
measure is not necessarily a good estimate of location Td of the 
ionization zone. 

Figure|6^ shows the value of -fd, as found from fitting the sig- 
nal in the frequencies, versus the value of Td, as determined from 
the location of the local minimum of Ti in the model. There is 
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Figure 6. (a) Plot of the fitted acoustic depth fd and (b) the correction term 
A(j, versus the acoustic depth as determined from the models and cor- 
responding to the local minimum in Ti (see Fig. 0. Filled symbols are 
for models using the CEFF equation of state, while crosses are for mod- 
els calculated using a simple Saha equation of state. The filled circles are 
for models having the same simple opacity (power law) but different theo- 
ries of convection, while the filled diamonds are models with the opacity at 
low temperatures from Kurucz. The values found for the solar data are also 
shown in both panels with 3<r error bars (dashed horizontal lines) due to the 
observational uncertainties. 



a difference of up to about 140 s between fd and Td, and one is 
not simply a function of the other. The difference between the two 
comes from a^, which measures the leading-order frequency de- 
pendence of the phase transition which the eigenfunctions undergo 
at the upper turning point. This will be strongly affected by the 
physics that change the structure of the surface, namely convection, 
EoS, the low temperature opacities, and the structure of the atmo- 
sphere. Consequently, we have to use some caution when taking the 
parameter fd from the fit to estimate the location of the ionization 
region. 

As an alternative we can consider one of the other parameters 
which also depends on the position of the ionization zone. This 
is Ad, given in Table [5] for all models and shown in Fig. |6j) as a 
function of the actual acoustic location of the ionization region. The 
value of Ad, defined in Eq. is not sensitive to the layers near 
the photosphere, as its value is determined exclusively by the sound 
speed at the ionization zone. However, the determination of this 
term is associated with a small correction in the amplitude, which 
makes it more sensitive to the observational errors when fitting the 
frequencies. 

Both panels in Fig.|6|show the solar values of fd and Ad with 
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Figure 7. Plot of the phase <f>o of the signal versus the envelope helium 
abundance Y for all models. The symbols are the same as in Fig.|f] The 
dotted lines illustrate the correlation between models with the same physics 
but different values of the surface helium abundance Y. The value found 
for the solar data is also shown with 3<r error bars (dashed horizontal lines) 
due to the observational uncertainties. 

3a uncertainties. The values of Ad indicate that all models calcu- 
lated with the CEFF equation of state give, even if marginally, a 
location for the ionization zone consistent with the Sun. 

Finally, the structure at the top of the envelope is also expected 
to be reflected in the value of 4>q. The value of this parameter for 
all models is represented in Fig. Q as a function of the envelope 
helium abundance. The larger difference is found when changing 
the EoS (about 0.06). But changes in the opacities also change <f>o 
by as much as 0.01, while the theory of convection changes this by 
about 0.01. It is interesting to confirm that the fitted value of (j>o 
is independent of the helium abundance, as one would expect from 
the analysis leading to the expression of the signal. Consequently 
4>o may allow the separation between the helium abundance and the 
physics relevant to the outer layers of the Sun because it is insensi- 
tive to Y whilst being indicative of some near-surface change that 
may be required in the physics. 

The solar value for 0o is also shown in Fig. Q Adjustments 
in the near surface layers seem to be necessary in order to produce 
models that have a value of 0o consistent with the Sun. Changes in 
the superadiabatic layer or in the surface opacities may be some of 
the options for reconciling the models with the solar data. 

4.2 The equation of state 

From the analysis of the results listed in Table|3| and as discussed in 
the previous section, the EoS is the most important factor in defin- 
ing the characteristics of the signal. In Fig. [8] we show the width 
parameter /3 as a function of Ad (a proxy for the location). Models 
that have the same EoS (CEFF) lie on a common locus in this dia- 
gram, as indicated by the dotted line. The position along this line of 
models all built with the CEFF varies according to changes in the 
convection or the surface opacities. Models Z\ and Z2, built with 
a different EoS, lie in a different region of the diagram. Thus we 
claim that, with the location of the ionization zone fixed, the width 
of the bump in Fi is mainly a function of the EoS, as expected. 

Consequently, after using 0o to test the surface physics, it is 
possible to combine the constraints provided by Ad and (3 to obtain 
a direct test on the EoS and the location of the ionization zone. 
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Figure 8. Plot of the estimated width of the bump, as given by j3, versus 
the value Aj providing an indication of the location of the ionization zone. 
The symbols are the same as in Fig.|6| The dashed line indicates a linear fit 
to the models with the same EoS. The values found for the solar data are 
also shown with 3<r error bars due to the observational uncertainties. 



Figure 9. Plot of the reference amplitude (see Eg. 1101 versus the envelope 
helium abundance Y . The symbols are the same as in Fig.|S| The dotted 
line illustrates the correlation between the helium abundance Y and the 
amplitude. 



Figure[5]also includes the parameters found for the solar data. 
These are marginally consistent with the expected behaviour found 
using models calculated with the CEFF equation of state. Other 
options for the EoS must be considered in an attempt to bring the 
models closer to the Sun. 

4.3 The helium abundance in the envelope 

From the discussion in the previous sections it follows that any de- 
termination of the helium abundance requires a careful tuning of the 
models to the correct structure of the envelope. Such a fine tuning 
can be performed based on the sensitivity of the eigenfrequencies 
to the behaviour of the adiabatic exponent in the region where he- 
lium undergoes its second ionization. We have found, as discussed 
above, that: 

• Ad provides a process to place the ionization zone in the 
model at the same acoustic depth as for the Sun - this corresponds 
to adapting mainly the surface layers of the model (atmosphere 
and/or convection) in order to place the ionization zones at the same 
acoustic location as measured in the Sun by the solar value of Ad; 

• j3 can then be used to adjust the EoS (or more likely to select 
it from a few candidates) to match the observed behaviour - this 
corresponds to verifying that the behaviour of f3 as a function of 
the location (Ad) in the models includes the observed solar values 
for these two parameters; 

• and finally, the parameters f<i and (j>o can be combined to 
adjust the surface physics in the model, in order to recover the 
observed solar values - this corresponds to adjusting convection 
(superadiabatic region, mainly), opacities (low temperature range), 
photosphere, etc, in a complementary way to the first point, until 
the solar values can be recovered with the models as both param- 
eters are strongly dependent on these aspects of the physics, but 
quite insensitive to the actual helium abundance. 

Consequently, we are left with one last parameter, connected 
with the shape of the bump through Sd, which is the amplitude of 
the signal ao, or A^. If the model has been adjusted to the observed 
data using the remaining parameters, then the amplitude will de- 
pend mainly on the helium abundance in the model, which can now 



be compared with the solar abundance. Such a relation provides 
a measurement of the helium abundance, which complements the 
boundary condition used in the evolution to fit the model to the 
present day Sun. 

Figure|5|illustrates how such a dependence of A<±, as defined 
in Eq. 1101 . could be constructed after the other aspects of the 
physics are adjusted. It is worth noting that, as expected from Fig.|5] 
the amplitude decreases with increasing Y, since the changes in 1\ 
due to changes on the hydrogen abundance dominate the variations 
of the bump. This regime for the inverse dependence of the ampli- 
tude of the signal on the abundance of helium is relevant for stars 
of low effective temperature. That follows from the overlapping of 
the three ionization zones (Hi, Hel, Hell ). For stars where these 
are fully separated in temperature it is expected that the amplitude 
will increase with the abundance of helium. 

As shown above (see Figs Q and |SJ the models used here are 
not fully consistent with the physics of the Sun and seem to be 
only marginally consistent regarding the equation of state that has 
been used. Consequently the amplitude A^, as found for the solar 
data, cannot yet be used as an indicator of the helium abundance 
in the solar envelope. A more adequate calibration of the surface 
layers in the models must be developed before an estimation for Y 
is inferred from this parameter. 

The simplified models we are using here to illustrate the ap- 
plicability of the method have been calculated with scaling a chem- 
ical profile determined without including diffusion and settling of 
helium. This is one of the aspects that needs to be considered in the 
models in order to move the parameters found for these closer to the 
solar values. With such a tuning, based on other seismic constraints 
and on the parameters of the signal discussed here, we have an in- 
dependent procedure to adjust our models to the Sun in this region 
near the surface, where the uncertainties in the physics dominate 
the structure of the models. 



5 CONCLUSION 

In this work we have developed a complementary method to con- 
strain the properties of the helium second ionization region near the 
surface of the Sun using high degree mode frequencies. The method 
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proposed here can independently test properties of this region, and 
provides a possible direct measurement of the helium abundance in 
the envelope. 

We have shown that some of the parameters characterizing the 
signature in the frequencies due to this region in the Sun are very 
sensitive to the EoS used in the calculation of the models, and so 
can be used to test and constrain the equation of state. Others of 
the parameters can also provide an important test on the physics 
affecting the surface regions of the models, namely convection and 
the low temperature opacities. By combining the diagnostic poten- 
tial of the five parameters determined from the data with very high 
precision the helium abundance can be effectively constrained. 

Here we were mainly concerned to establish the method and 
demonstrate how it can be used to study the Hell ionization zone 
in the Sun, and the physics that affect the structure of the Sun in 
that region. In spite of having used simplified models to represent 
the Sun we have illustrate the sensitivity of each parameter to the 
physics, establishing the approach that can be followed when ade- 
quate up-to-date evolved solar models are used. Besides the phys- 
ical ingredients addressed here, aspects like diffusion and settling 
and improved opacities have to be implemented in order to provide 
a physically consistent value of the helium abundance. A calibra- 
tion of the actual solar helium abundance using models with the 
best up-to-date physics will be the subject of the second paper in 
this series. 
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APPENDIX A: A VARIATIONAL PRINCIPLE FOR THE 
Hen IONIZATION ZONE 

W e consider here a variational principle, following the procedure 
by Mo nteirol ll99d) . for describing how the modes are affected by 
the presence of the region of the second ionization of the helium. 
We start by using a variational principle, for small changes of the 
eigenfrequenci es (ui) due to small chan ges of the structure. It can 
be written (see lChristensen-Dalsgaard et al Jl 19951 and references 
therein) in the form 

ouj = — 



h 



with 



Jl ~ ~ r t E 2 



Here, r t is the acoustic size of the Sun, and 
d8B Q \ „i . , „ dE 2 



SI 



J 



A 2 E 2 , 
dr 2 



(Al) 



dr , (A2) 



where E r is the normalized radial component of the eigenfunction 
(with constant amplitude Eo). The acoustic depth r is defined in 
Eq. J2}. 

From asymptotic analysis (see MCDT) we also have that well 
inside the turning points and for moderate degree modes, 



<TE r 
dr 2 



-UJ (1-A) E r 



E r ~ Eo COS 



(l-A) 1/2 dr + 



(A3) 



(A4) 



The changes in the structure relative to the reference ('smooth') 
model are described with the functions SBi, as given by 
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(A p -A c ) 2 
4A 



(1-A) 2 J (TiP) 
1-2A 



A„ + 
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+ 



(1-A) 2 
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4A 



-A 



8B 2 
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+ 



o 1-3A/2 
(1-A) 2 
1-2A 



+ 



5p 

1-A A p -Z 
2(1-A) 2 2A 9 



(A6) 



(riP) 



+ ■ 



A 



(1-A) 2 ' (1-A) 2 
«(TiP) , 1 



A p -A c 
2A D 



A 



2 1-A (fiP) 2 (1-A) 2 



5p 
P 



(A7) 



(A8) 



where r, p, c and 5 are distance from the centre, density, adiabatic 
sound speed and gravitational acceleration, respectively. We have 
also introduced the following quantities 

l(l+l)c 2 



A 

where I is the mode degree, and 

A P 



g d g 
—2~ -5- log — 

oj 2 c dr \ pc 



A c 



A Q 



uj 2 c dr 
(-) 2 



(A9) 

(A10) 

(All) 
(A12) 



These are all first order quantities, compared to unity, because well 
inside the resonance cavity of the modes the local wavelength is 
significantly smaller than the scale of variations of the equilibrium 
quantities. 

In order to use the expression for 81 from Eq, IA2l it is nec- 
essary to replace the term in (dSBo /dr) by integrating by parts to 
obtain for 81; 



81 



SB 1 E 2 r + (SB 2 +SB ) ^ + 8B 3 ^f 
dr dr 2 



dr. (A 13) 



The integration is done only for the region of the ionization zone, 
starting at r a and ending at Tb. Because we are restricting our anal- 
ysis to localized variations, it is also assumed that the model dif- 
ferences are zero everywhere else. This does not affect our result 
since we will only take those changes in the frequencies that are 
not affected by model differences spreading over regions of size of 
the order of (or larger than) the local wavelength of the modes. 
We recall, from asymptotic analysis, that 

p2 



AEt 



\ cos(A) 
" El 



M1~A)^ sin(A), 



(A14) 



^ 4c 2 (l-A) cos(A) . 
The argument of the trigonometric functions is 

A(r) = 2 lu ( (1-A) 1/2 dr + <j> . (A15) 
. Jo 

After replacing these expressions in the equation for SI, we find 



S 2 E 2 



81 



SB! 



■± -4(1-A)<5£ 3 



cos A 



_ 2(1 _ A) V 2 ^±^ sinA l dr . ( A16) 



This expression gives the variational principle for perturbations in 
the frequencies due to small changes in the structure, as described 
by SBi. 

The next step is to establish what is the effect on the structure 
of the ionization zone for helium, relative to a model where such 
a localized effect is not present. In particular, we need to estimate 
how Ti, P and p are changed from being slowly varying functions 
of depth to the actual values they have when the second ionization 
of helium occurs. The difference will produce the 8(TiP) and Sp 
responsible for changing the frequencies, as given in Eos IA5IA8I . 
That will allow us to calculate an expression for the characteristic 
signal we want to isolate in the frequencies. 

In order to find an expression for the signal we shall first con- 
sider that the changes are dominated by Ti. In doing so, we adopt 
here a different approach from lMonteirol fl996). who consider that 
the dominant contribution could be isolated in the derivative of the 
sound speed. We do so because the effect of the ionization is better 
represented as a 'bump' in Ti (see FigsQand|2j, extending over a 
localized region of the Sun. Therefore we retain the terms for SVi, 
and neglect, as a first approximation, the contributions from Sp and 
SP. In doing so we assume that the changes in the sound speed are 
mainly due to the changes in the adiabatic exponent. 

Now, relating SI to the change in the eigenvalue Suj (and using 
Eq.IAII it follows that 



[5u 



\SIl 



iri 



f Tb 8Y 
— / (/ c cosA + / s sinA) — - dr , (A17) 
2r t / Ti 



where / s and f c are functions obtained from adding the coefficients 
of STi in the expressions of SBo, 8B1, 8B2 and 8B3 (see Eq. lA16l 
and Eas lA5IA8> . 

At this point we introduce an approximate description of the 
effect of the second ionization of helium on the adiabatic expo- 
nent. As represented in Fig.|2p, we adopt a prescription where the 
'bump' is approximately described by its half width and height 
Sd = (STi/ri) Td , with the maximum located at ra. This corre- 
sponds to considering the following approximating simple expres- 
sion for 8T1: 



STi 



= <5d < 



1- 







r— r d 


T-T d 







; Td— (1— a)(3 < r < r d 

; r d < r < r d +{l+a)l3 
; elsewhere. 



(A18) 



The region of the ionization zone starts at r a =rd — (1— a)f3 and 
finishes for Tb=Td+(l+a)f3, giving that Tb— r a =2/3 is the width. 
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The parameter a represents the asymmetry of the bump, and for a 
first order analysis it does not affect the result. 

We further consider that the functions f 3 and f c are slowly 
varying functions of the structure when compared with the size 
of the ionization zone (~2/3), and so their derivatives can be ig- 
nored in the integration. Using this approximation we may integrate 
Eq. IA17i finding that 



2r t r | w/3(l-A)V2 

cosA d + f s sinA d ) . (A 19) 



(fc 



All quantities are now evaluated at T=Td- 

Taking the dominant contributions (in terms of powers of ui 
and derivatives of the reference structure - see CDMT for details) 
of the functions f c and f s (Eq, IA16t . we can finally write the signal 
as being 

. 35 d 1-2A/3 sin 2 [^(l-A) 1 ' 2 ] 
[Hri ~2^^=A Z^9(!=A) cosA - ( A20 > 

This is the expression that describes the 'additional' contribu- 
tion to the frequencies of oscillation u> n i if the region of the second 
ionization of helium is present. By assuming that we have 

Unl = [tOnl] smooth + [5w n ;]r! , (A21) 

it is now possible to try removing the smooth component 
[w„z]amooth, by adjusting the frequencies to the expression we have 
found for the 'periodic' component [(Java]]?!. In doing so the pa- 
rameters describing the structure of the Sun at the location ta are 
determined. 



